function f = getS(S,N,share,delta,alpha,psi)

A = getAm(N,S,alpha,psi,delta);
f = abs(A.*share - psi.^(S./(1-alpha)));

end
